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Abstract. We consider a system of A'' classical particles, interacting via a smooth, short- 
range potential, in a weak-coupling regime. This means that N tends to infinity when the 
interaction is suitably rescaled. The j-particle marginals, which obey to the usual BBGKY 
hierarchy, are decomposed into two contributions: one small but strongly oscillating, the 
other hopefully smooth. Eliminating the first, we arrive to establish the dynamical problem 
in term of a new hierarchy (for the smooth part) involving a memory term. We show that 
' ^ '■ the first order correction to the free flow converges, as N ^ oo, to the corresponding term 

^ ^ associated to the Landau equation. We also show the related propagation of chaos. 

1. Introduction 

Lev Landau in 1936 proposed a kinetic equation, usually called Fokker-Planck-Landau 
equation (simply Landau equation in the sequel) which is a diffusion with friction in velocity, 
suitable to describe the behavior of a weakly interacting gas, in particular a Coulomb gas in 
a regime where the grazing collisions are dominant. 

Roughly speaking the Landau's argument was to take the Boltzmann equation with Coulomb 
cross-section and (cutting-off short and long distances) apply the Taylor expansion to the col- 
lision operator. The result is a degenerate elliptic operator acting on the velocity space (see 
[18] and the original publication of Landau [17]). The full Taylor expansion of the Boltzmann 
collision integral for arbitrary intermolecular forces was studied in [6] and a formal general- 
ization of Landau collision integral to arbitrary scattering cross-section was proposed there. 
A more precise asymptotics in the Coulomb case was also studied in [8]. 

The Landau equation for the one particle distribution f{x,v,t), where x G M^, f G and 
t € M+ denote position, velocity and time respectively, reads as 

(1) idt + v-V,)f = QLif,f) 
with the collision operator Ql given by: 

(2) QUf, f){v) = I dviV, [a{v - vi) (V, - f{v)f{vi)] . 



Here x plays the role of a parameter and hence its dependence is omitted. Moreover the 
matrix a{w) has the form 
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(3) 



a{w) 
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where ^4 > is a suitable constant. 

Note that the Landau equation possesses all the properties known for the Boltzmann equa- 
tion, namely the mass, momentum and energy conservation and the H-thcorcm. Actually the 
homogeneous Landau equation can be rigorously derived in the grazing collision limit of the 
homogeneous Boltzmann equation by a suitable rescaling of the cross-section. 

In particular, in [1] the authors show that, under suitable assumptions on the cross-section, 
the diffusion Landau equation (1) can indeed be derived. The diffusion operator is the form 
(2) but with a matrix a replaced by 



with a a smooth function. Next in [13] and [23] steps forward were performed to arrive to 
cover the case ~ for small 1^1, with u <\. 

The case of the matrix (3) was treated in [24]. It is worth to underline that the initial value 
problem for the homogeneous Landau equation is strongly simplified for the case ad^yj) ~ j^^j 
with v < 1 (see [9] and [10]), while for the matrix (3) we have a weak existence theorem 
obtained by compactness arguments based on the entropy production control [24]. Moreover, 
for the inhomogeneous case, we have existence and uniqueness of strong solutions for data 
sufficiently close to a Maxwellian [14]. This is the only existence and uniqueness result we 
are aware. 

A natural question is to see whether the Landau equation can be directly derived, under 
a suitable scaling limit, from a particle system as it is the case of the Boltzmann equation. 
In fact one can sec ([2], sec also [22] and [21]), at a formal level, that the Landau equation 
is expected to be valid for a weakly interacting dense gas. The precise statement and scaling 
(called weak-coupling limit) will be presented and discussed in the next Section. The formal 
analysis gives indeed the Landau equation (1) with matrix (3). The two-body interaction 
potential (f) is assumed smooth, spherically symmetric, and the constant A is given by: 



where 4>i\k\) = J dx (j){\x\)e~^'''^ . 

Note that we find the Landau equation with matrix (3) , which is not related to the Coulomb 
potential, but arises even though the potential is smooth and short-range. This fact was first 
established by N.N. Bogolyubov in 1946 [20]. 

In the present paper we want to start the rigorous analysis of the weak-coupling limit for 
an Hamiltonian particle system. Our result is very preliminary. We first decompose the j- 
particle marginals into two terms, one hopefully smooth and the other strongly oscillating. 
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but small. Eliminating this last term from the equations (with a procedure similar to that 
proposed by Zwanzig [26]) we find an equation with memory, which we can handled up to the 
first order in time. We show that this contribution agrees with the corresponding one arising 
from the Landau equation. Roughly speaking we present a rigorous derivation of the Landau 

equation at time zero. 

It is well known that the situation for the Boltzmann equation is better, namely we are able 
to derive such a kinetic equation for a short time [16] (see also [7] for additional comments 
and results) in the low-density (or Boltmann-Grad) limit. 

Note that the linear case, namely a single particle in a random potential under the weak- 
coupling limit, is well understood, see [11] and references quoted therein. 

Our analysis deals with the nonlinear problem but our techniques could apply as well to the 
linear case. Wc think that, while we can easily obtain the same consistency result presented 
here, it seems very difficult to go further. In [11] and related references, it is crucial the use of 
probabilistic tools which seems more efficient compared with the hierarchical approaches. In 
contrast it is very difficult to implement the ideas working for the linear case to the present 
problem. 

Finally we want to mention that the same problem of characterizing the weak-coupling 
limit of particle systems, arises also in a quantum mechanical context. In this case the 
quantity which we are interested in is the Wigner transform [25] which is a way to describe a 
quantum state as a function in the classical phase space. In contrast with the classical case, 
we expect that the Wigner transform approaches, in the weak-coupling limit, the solution of 
a suitable Boltzmann equation, with a corrections due to the statistics, whenever taken in 
explicit consideration. We quote [15], [3], [12], [4], [5] for the few results in this direction and 
[19] and references quoted therein, for the Boltzmann description of wave dynamics in the 
weak-coupling limit. 

2. Weak-coupling limit for classical systems 

We consider a classical system of N identical particles of unit mass in the whole space. 
Positions and velocities are denoted by the vectors Qn = {qi . . . qn} and Vn = {vi . . . vn} 
respectively. The particles interact via a spherically symmetric, smooth potential of finite 
range <^ : — >■ R, namely (f){x) = if |x| > r for some positive r. In the following we assume 
units for which r = 1. 

The Newton equations read as: 

(5) -^qi = Vi T"^^ = X] 

^ j=l...N: 

Here F = — V(/) denotes the intcrparticlc (conservative) force, and r is the time. 

Let e > r be a small parameter denoting the ratio between the macroscopic and microscopic 
space-time unities. 
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We axe interested in a situation where the number of particles N is very large and the 
interaction strength quite moderate. The system has a unitary density so that we assume 
N = e~^. In addition we look for a reduced or macroscopic description of the system. Namely 
if q and r refer to the system seen in a microscopic scale, we rescale eq.n (5) in terms of the 
macroscopic variables 

X = eq t = ET 

whenever the physical variables of interest are varying on such scales and are almost constant 
on the microscopic scales. 

Remembering that we want to describe weakly interacting systems, we also rescale the 
potential according to: 

(6) (t> ^ V^(/>, 

so that system (5), in terms of the variables, becomes: 

^ j=l...N: ^ j=l...N: 

Note that the velocities are automatically unsealed. 

A statistical description of the above system passes through the introduction of a probabil- 
ity distribution on the phase space of the system. Let = W^{Xn, V^v) be a symmetric (in 
the exchange of variables) probability distribution. Here {Xn, Vn) denote the set of positions 
and velocities: 

Xn = {xi...xn} Vn = {vi...vn}, XieW^,VieR^. 
Then from eq.ns (7) we obtain the following Liouville equation 

^ 1 

(8) (dt + V,JT^^(X^, Vn) = [TlfW'') {Xn,Vn). 

i=l ^ 

Here we have introduced the operator 

(9) {T%W''){XN.Vn)= J2 {n,iW''){XN,VN), 

0<k<e<N 

with 

(10) T^^,W^ = V<^(^^) • (V., - V.,)W^. 

To investigate the limit e — > it is convenient to introduce the BBKGY hierarchy for the j- 
particle distributions defined as 

(11) fj^i^j: ^j) = J dxj+i ...J dxN J dvj+i ...J dvN 

W^{Xj,Xj+i . . . xn; Vj,Vj+i . . . vn) 
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for j = 1. . . . , AT - 1. Obviously we set = W^. Note that BBGKY stands for Bogolyubov, 
Born, Green, Kirkwood and Yvon, the names of physicists who introduced independently this 

system of equations (see e.g. [2]). 

Such a hierarchy is obtained by means of a partial integration of the Liouville equation (8) 
and standard manipulations. The result is (for 1 < j < N): 



(12) 

for 1 < i < AT. 
We set 



/f = 0, for j>N, and = W^^. 



The operator C?^^ is defined as: 

j 
fe=i 

and 

(14) Clj+Jj+i{xi ...Xj;vi... Vj) = 

- J dxj+i J dvj+iF ^^^__3±ij ■Vv^fj+i{xi,X2,. . . ,Xj+i;vi,. . . ,Vj+i). 

Cfe j+i describes the interaction of particle k, belonging to the j-particle subsystem, with a 
particle outside the subsystem, conventionally denoted by the number j + 1 (this numbering 

uses the fact that all the particles are identical). 

We finally fix the initial value {fj}f^i of the solution {/j^(t)}jLi assuming that {/j'jjLi 
is factorized, that is, for all j = 1, ... A'' 

(15) /° = /o^^ 

where /o is a given one-particle distribution function. This means that the state of any pair of 
particles is statistically uncorrelated at time zero. Of course such a statistical independence 
is destroyed at time t > because dynamics creates correlations and eq.n (12) shows that the 
time evolution of ff' is determined by the knowledge of which turns out to be dependent 
on and so on. However, since the interaction between two given particles is going to 
vanish in the limit e — )■ 0, we can hope that such statistical independence is recovered in the 
same limit. Therefore we expect that when e — > the one-particle distribution function 
converges to the solution of a suitable nonlinear kinetic equation /, which we are going to 
investigate. 

If we expand /j^(t) as a perturbation of the free flow S{t) defined as 

(16) iSit)f,)iX„V,) = f,{X, - V,t, Vj), 
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we find 

(17) //(*) =s{t)f^ + s{t - h)q^j^'i,{h)dh+ 



^ f S{t-h)T^f^{h)dh 
Jo 



We now try to keep information on the limit behavior of fj{t). Assuming for the moment 
that the time evolved j-particle distributions fj^{t) are smooth (in the sense that the first 
and second derivatives are uniformly bounded in e), then 

(18) q^jll,{Xj;Vj;h) = 

-e^^ J dr J dvj+iF{r) ■V^^fj+i{Xj,Xk- er;Vj,Vj+i,ti). 

k=l 

Because of the identity 

(19) J drF{r) = 0, 
we find that 

(20) q+J^_,,{X,;Vf,h) = 0{e^) 
provided that D^fj^^ is uniformly bounded. Since 

N-j ^, 7, 

we see that the second term in the right hand side of (17) does not give any contribution in 
the limit. 
Moreover 



(21) f S{t-ti)T^f^{h)dh 

Jo 



where / is a smooth function. We note that the time integral in (21) is 0(e) because F ^ 
only for times in an interval of length 0{e). Therefore /j^ cannot be smooth since we expect 
a nontrivial limit. 

In order to look for a (nontrivial) kinetic equation, we can conjecture that 
(22) /f = 5f + 7f 

where g^^ is the main part of /j^ and is smooth, while 7^ is small, but strongly oscillating. 
We operate this decomposition according to the following equations which define and jj^: 

3 



(23) {dt V. j^f = ^q+i5f+i + ^q+i7f+i 
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(24) {dt + p^v,. V. J7f = ^T^'lf + ^T^g^, 
with initial data 

(25) g^{Xj,Vj,0) = f^{Xj,Vj), jf{Xj,Vj) = 0. 

Note that 7f = since Tf = 0. 

The remarkable fact of this decomposition is that 7 can be eliminated. Indeed, let 

iX,it),Vjit)) = i{xiit)...xj{t),vi{t)...vjit)}) 

be the solution of the j'-particle flow (in macro variables) 



d d 
26 — Xj = Vi —V. 

^ ' dt dt 

-3 



1 ^ v^f^t^), 



with initial datum (Xj, V^) = {\x\ . . . Xj,vi . . . Vj}). Denote by Uj{t) the operator 

(27) Uj{t)f{Xj,Vj) = exp{t(-^^;, • V., + -^Tj)}f{Xj,Vj) = f{Xj{-t),Vj{-t)), 

then eq.n (24) can be solved: 

(28) 7f (*) = f dsU,{s)^T,gf{t - s). 

Jo 

Explicitly 
(29) 

7f(x,,y„t) = -i= Td. J2 F ( ^-^-^^ 7^^"^^ )-[(v. -v.jgf](x,(-.),y,(-.);t-.). 

l<i<k<j ^ ^ ^ 

Inserting (28) in (23) we finally arrive to a closed hierarchy for {ffj^j^Li- Obviously we pay 
the price of a memory term given by the time integral in (28) or in (29). 
We write the hierarchy in integral form. Then 

(30) g^{t) =S{t)f^ + ^ f S{t - T)q^,gf^,{T)dT 

Jo 

+ ^ f S{t-T)q^,^f^,{T)dT 

Jo 

=S{t)f^ + ^ 1^ S{t - T)q+,gf+,{T)dT 

+ f dr r daS{t - T)q+,U,+i{T - a)T,+,g^^,{a). 

£ Jo Jo 



8 



A. V. BOBLYLEV, M. PULVIRENTI, AND C. SAFFIRIO 



Remark 1. Why do we expect that strongly oscillates? Let us try to control the first 
derivatives ofh{Xj, Vj,t) = Uj{t)ho{Xj, Vj) = ho{Xj{—t),Vj{—t)) for a given smooth function 
Hq . Then 

"t^te^ ^^^ ^^^^^^^ ^^^ ^^^^j 



and analogous formula for . Here we are using Greek indices for the components of Xi 
3 quanta 

( changing —t-^t) 



and Vi. To estimate quantities like ^^gja , —§^cr^, ^^dv" ™^ ^^'^ (^^) '^^^ fi''^^ 



ddxjjt) _dvl{t)^ 

dt dx? dx? ' 



(31) 



d Hit) _ 1 V dF^ xk{t)-xr{t) ( dxlit) dxUt) 
^ ' dt dxf £3/2 2^ dxl^ £ '\ dxf dxf 

r^k 

Integrating eq.ns (31) and (32) in time, we arrive, by using the Gronwall lemma, to 



dvlit) 



dxf 



< C exp 



Ctc 

.3/2 



where Tc is the scattering time, namely the time interval for which \xk{t) — Xr(t)\ < e . Now, 
even though Tc = 0{e) (neglecting small relative velocities), it seems difficult to get something 
better than a bound like exp{-^). 

In conclusion we expect that the first derivatives of h{t) are 0(exp(-^)). Looking at eq.n 
(28) we expect for 7 the same behavior. In contrast, the action of the operator Cj is regular- 
izing (althoug we are not able to prove this) so that we expect g to be smooth. 

On the other hand is also expected to be small, in some sense. Indeed by taking the 
scalar product of (28) by a smooth function u, we find 



(u,7f (t))| < ^ f ds\\Uj{-s)u\\L^\\T^g^{t-s)\\Li 
Jo 



< £ 



5/2 



|(V„i -Vv2)9j {xi,xi + £r,X3...,Vj-t- s)\. 
Therefore this term is vanishing provided that g^ is sufficiently smooth (uniformly in e). 

A rigorous analysis of the limit N ^ 00, e = A^~(^/3) seems to be very difficult. We expect 
that, in this limit, both fj^{t) and gf{t) would converge to f{t)®^, where / solves the Landau 
equation stated in Introduction. We cannot prove it, but a first step in this direction is made 
in the following Sections. 
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3. CONSISTENCY 
We consider eq.n (30) written in symbolic form as 
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(33) 5, = S(t)/0 + ^,+i5,+i, 

where all upper indices N are omitted for brevity. To solve these equations one can use the 
obvious iterative scheme 

9^ = 5(i)/°, 5^ = S(i)/° + Aj+ig^%, n = 0, 1, . . . 

Our goal in this section is to prove that the equation for g\^\t) = gi{t) is consistent with 
the Landau equation. Thus we replace (30) by its first approximation: 

(34) g^it) =S{t)f^ + ^ f S{t - T)q^,S{T)f^^,dT 

Jo 



(35) =S{t)f^ + ^ f S{t - T)q+,S{T)gf^,dT 



+ !L^ f dr rdaS{t~T)q^,Uj+,{T-a)Tj+,S{a)f^^,. 
^ Jo Jo 



Here we set 



(36) ^f{Xj, Vj,r) ~ r daUjir - a)TjS{a)f^^, 

Jo 

(37) =-^fds V p^ -i{-s)-x,{-s) ^_ 

• [(V,, - V,J5(r-.)/0](X,(-.),y,(-s)). 



We note that can be explicitly computed. 



Lemma 1. We have 



^^{Xj,Vj,t) = {Uj{t)f^ - S{t)ff){Xj,Vj) 
(38) = ff{Xj{-t), Vj{-t)) - f^{Xj - Vjt, Vj) . 
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Proof. Let = — Y^- Vi ■ V^. be the free flow generator. 
Then we compute 

(39) U,{t)f^ - S{t)ff = 1^ ds-^^[U,{s)S{t - .)]/° 

= fds[Uj{s)iLo + -^T,)S{t - 

- / ds[Uj{s)LoSit - s)]f^ 
Jo 



For convenience of the reader we make exphcit eq.n (34) in the case j = 1 

iV - 1 /•* 
e Jo 

N-1 



(40) g^it) =S{t)fo + / S{t - T)CIS{T)f^dT 

^ J^S{t-T)Clj^{T)dT 



+ 

where, by Lemma 1, 



1 /"^ 

(41) ^^{xi,vi,X2,V2,t) = -^ dsF{ 

V £ Jo 



1 r , ,Xl{-s) - X2{-S) 



£ 

■ [(V„, - V,,)S{t - s)f^]{X2{-s),V2{-s)) 
= [/2°(X2(-r), V2(-t)) - /2°(X2 - V2r, V2)] . 
The first result of the present paper is summarized in the following Theorem. 
Theorem 1. Suppose /o G Co(M^ x R^) he the initial probability density satisfying: 

(42) |D7o(x,^;)| <Ce-^l''l' for r = 0,1,2 

where is any derivative of order r and b> 0. Assuming also that cj) G C'^{E?) and (j){x) = 
if\x\ > 1. // (15) holds for j = 1,2, then 

(43) lim 5f (i) = S{t)fo + f dTS{t - T)QL{S{T)fo, S{T)fo), 
£->o Jo 

(44) lim7f(i) = 0, 

£->0 

where Ne^ = 1 and the above limits are considered in D' . 

Proof. Let u G D{M.^ x R^) be a test function. Prom now on we will denote by {hj,kj) = 
/ dXj J dVj hj{Xj,Vj)kj{Xj,Vj) the inner product. Then 

(45) {u,g^{t)) = iu,Sit)fo) + ^^ / iu,Sit-T)CIS{T)fl)dT,+ / T,(r)dr, 

Jo Jo 
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where 



(46) 7e{r) = - ^ ^ J J J dvi J d'U2(V^i5(r - ui))- 



Xl-X2\ _7V/ N 
72 {X1,X2,V1,V2,T). 



We have already seen that the second term in the right hand side of (45) is vanishing. There- 
fore we have to evaluate the last term, namely Jq dr'j£(T). We split the term 7^{t) into two 
terms 

(47) T, = T|+T> 
where 

(48) J-^(r)>=-^^^ / dxi [ dx2 I dvi [ dv2{Vv,S{T - t)u{xi,vi)y 

J J J J\w\>ae^/* 

/xi — X'2\ ^TV/ \ 
■Fi 1 72 [Xi,X2,Vi,V2,T) 

where w = vi—V2is the relative velocity and a is a number to be fixed later on. 7f is defined 
accordingly. 

The reason of this decomposition will be clear later on. For the moment we show that 7f is 
negligible. 

Lemma 2. 

(49) 7f = 0{e'/') . 

Proof. By Lemma 1 we have that 7^ is uniformly bounded. Moreover by the change of 
variables 

X2 = xi — sr 

we get 

(50) \7f\ <C{N -l)e^^ [ dxi [ dvi\Vy,S{T - t)u{xi,vi)\ [ dr\F{r)\ [ dw 

J J J i|«)|<0£V4 



To evaluate we use (41) to write it as 
(51) 7>{t)={N -l)e^ I dxi f dr f dvi f dv2 

J J J J|w)|>a£l/4 

1 /"^ , ^ , , ^ / Xli — s) — X2{ — S)\ , 

- dsFa[r)F0\ 1 [he{xi,X2,Vl,V2,T,s)\a,p, 
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where X2 = xi — er and is the matrix 

{he)a,fi = -C^v^Sir - t)u{xi,v,))a[{Vy, - V^,)S{t - s)f^UX2{-s),V2{-s)), a, /? = 1,2,3. 

The summation over repeated Greek indices is assumed here and below. 

Here the flow X2{t) = {xi{t),X2{t)) has initial conditions {xi,xi — er). Scahng times we also 

find 



dV2 

'|i«|>a£l/4 



(52) 7>{t) ={N-l)e^ [ dxi [ dr [ dvi [ 

J J J J\w\>ae'^/ 

1 / ( x\{—es) — X2{—ss)\ , 

y dsFa.{r)F^ 1^ 1 V^z\X\,X2.,Vx,V2,T,£s)\o.,^ . 

Let us introduce the function h which is the formal limit of /ig, namely 

(53) ho,,fi(xx,vx,V2,T) = -i?a(a;i,vi,r)[(V^i - V^2)5(r)/2(xi, 2:1,^1,^2)]^, 
where 

(54) R(x\,v\,t) =Vvj^S{t - t)u{xi,vi). 
We spht 7^ into two terms 

7f = 7i + 

where 

dV2 



(55) 



7>{t) = {N-l)e^ [ dxi [ dr [ dvi [ 

J J J J|«,|>ael/4 



f ^ , ^ , , „ / xi(—es) — X2(—ss)\ , 

dsF^{r)Ffs i-^ '-j^ n ha,^{xi,vi,V2,T) 



and 

dV2 



(56) 



7>{t) = {N- l)e^ [ dxi f dr f dvi [ 

J J J J|t«|>ael/4 



We shall show that (r) is vanishing while 7^ (r) has the right behavior. In the evaluation 
of (r) we note that h does not depend on s so that we have to evaluate the integral 

(57) IJ dsF j^^iizffl^^^) = i dsF ^M-s)-X2{-s) 

Indeed the integral (57) can be bounded when the interaction time of the two-particle system 
is 0(e) and this is true only if the relative velocity is not too small (see Lemma 3 below). 
This explains why we did the decomposition (47) . 
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Lemma 3. Setting w = vi — v-2, suppose that 

(58) > ae^/^ 
where a = 4y^||F||x,oo . Then, defining for any real number s 

(59) A, = {s\\xi{s)-X2is)\<£}, 
we have 

(60) meas{A^) < --. 



\w\ 



Moreover, for i = 1,2; 

(61) \vi{es) - Vi\ < C^. 

Proof Assuming first that s > 0, we pass in the coordinate system around the center of mass 
(at the origin) and denote by ^{t) = xi(t) — X2{t). Let w = vi — V2 he the relative velocity and 
Wx its horizontal component. We assume that at time zero the particles are in the interaction 
disk (more precisely, they enter in the interaction disk at time s = 0) and fix the axis in such 
a way that w is horizontal and its x- component is positive, namely Wx = \w\. Let t be the 
first time for which 

Wx{t) < —. 



By the equation of motion 
(62) 

we infer 
from which 

(63) t> ^ ' ' 



Wx{t) = \w\+ J -^^x ^ 



— >\w\ p F Loct 

2 s/e 



4||F||z,oo 

In the time interval [0,t\ we have Wx > ^^'i horizontal displacement is (under assump- 
tion (58)) larger than 

(64) ^*'^2£, 

since the diameter 2£ is a maximal path inside the sphere, independent of the initial point. 
This implies that, when \^{t)\ < e, then \w{t)\ > |u;(0)|/2 and hence 

4e 

(65) meas(A£) < - — -. 

\w\ 



Moreover 

/•es 

lU ( \ =71-1-4- / 



(66) v^{ss) = v,+ l ^f( "^^"^""^^"^ ) da 
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(67) 
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\vi{es) — vi\ < 

\w\ 



The case s < reduces to the case s > by changing the initial velocities to Vi{0) = — ■Uj for 
i = 1,2. This completes the proof. ■ 



Note that 



(68) 



xi{—es) — X2{—es) 



1 f~^^ 

= r — ws + - / da[{vi{a) — vi) — {v2{(t) — V2)] 
^ Jo 



xi{—es) — X2{—£s) 



(r — ws) 



thus, by Lemma 3, 
(69) 

The integral (55) reads 

(70) (r) = {N - l)e^ J dxi j dr j dvi J 



<CsA 



\w\ 



dV2 



'|«)|>aeV4 

dsFa{r)Fjs{r - ws)ha,i3{xi,vi,V2,T) + E 
where the error term E is given by 



(71) E = {N- l)e^ [ dxi f dr f dvi [ 

J J J J\w\>a6^ 



dV2 



Jo 



Fa 



|tu|>a£V4 

xi{—es) — X2{—ss 



— Fp{r — ws) 



ha,p{xi,Vl,V2,T) 



It is clear from the proof of Lemma 3 that \xi{—es) — X2{—£s)\ > e if s > 4/|u;| (see (64)). 
On the other hand, \r — ws\ > 1 if s > 2/|u;|, provided \r\ < 1. Hence, 



(72) \E\ <Cy/e [ dxi [ dvi [ '^'"27—7 /' ' sds\h{xi,vi,V2,T)\ 

J J J\w\>ae^/^ H Jo 

<C^ [ dxi I dvi [ dv2T^\h(xi,vi,V2,r)\ 

J J J\w\>ae^/^ H 



<C^/e\ \oge\. 
In the last step we estimated 

(73) \h{xi,Vi,V2,T) <C|/o(a::i - V2T,V2){'Vvi -TVa:i)fo{xi -ViT,Vi)\ 

+ \ fo{xi -ViT,Vi){Vv2 -rVx2)fo{xi -V2T,V2)\ 
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Lemma 4. For all w j^O, 

/r+r/e f f+^h 

dr I dsFa{r)Fp{r — ws) = - lim dr dsFa{r)F^{r — ws) = a{w)a 

Jo 2e-^0j J-r/e 

where 

(75) aiw)a,(S = J^i^a,(5 - 

and 

1 f°° 

(76) ^ = ^1 dpp'cl^^p), 







with 0(1 /c|) = /jj3 0(r)e- 



-ik-r 



Proof. The first identity in (74) is due to tlie symmetry F{r) = F{—r). Then we compute 
tlie left liand side of (74) taking the Fourier transform and passing in spherical coordinates. 
The result is 

(77) A dk5{k-w)k®k = a{w). 



Finally by the use of the dominated convergence theorem we can establish 

(78) hm / dT7>{T) = f drSit - T)QL{S{T)fo, 5(t)/o) 
£->o J Jo 

in D'. To conclude the proof it remains to show that 

(79) lim / dT7^{T) = 0. 

We first evaluate 
(80) 

{he - h)a,0{xi,r,Vi,V2, T, £s) = i?„ (xi , -Ui , T ) 

Note that 

V,5(t)/(x, v) = 5(t)(V^ - rV^)/(x, v). 
Omitting irrelevant variables we observe that 

{he - h)a,f} = Ra{^f}{-es) - $^(0)) 

where $(a) = [(V^^ - V^J5(r + u)f^]{X2{a), V2{a)). 
Hence ^ 

\he-h\ < \R\ I d(j\^{a)\. 

J —es 

It is easy to see that ^(cr) is a linear combination of various second derivatives of , multiplied 
by w{a) = -^F ^EiMzE^M^ ^ pi^s two terms proportional to first derivatives with respect to 
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X. All the derivatives are computed at the point [X2{a) — {t + a)V2{cr),V2{a)]. Hence, under 
the assumptions of Theorem 1, we obtain 

\he -h\< C\R\^ / daexp{-6(|i;i(-a)|2 + \v2{-a)\^)}. 
v£ Jo 

Since \xi{—es) — a;2(— £s)| > e if s > 4/|'u;|, the integral over ds in (56) can be estimated by 



(81) 



, ^ fxi(-es)-X2(-es)\,^ C\R\ /"^/l"'! , , 

dsFp i-^ L—^ '-j {he - h)a,p < ds da^ia) 



C\R\^ 



i/\w\ 



< —. — i— ^ / daip(ea) 

m Jo 



where ip{a) = exp{— 6(|vi(— (t)P + |t'2(— c)!^)}. 
Then by energy conservation 

\viit)f + \v2it)f + 2^cl> j^ ^i(^)-^2(0 ^ ^ ^^^^^ 

and therefore 

iPies) < exp{-b{\vi\^ + \v2\^ - A^UW^)}. 
Hence, we obtain the following estimate of the integral (56): 

\7>{T)\<C^/^ J dxi j dr j dV2\vi-V2\-''\R{xi,vi,T)\\F{r)\ exp{-6(|i;i|2 + |^;2|')}, 

where R{x\,vi,t) is given in (54). It is clear that R{xi,v\,t) = if > Ri, where Ri 
depends only on r. Therefore 

\7^ir)\ < C^e J dvi J dv2\vi - V2\-'^ expi-bUvil"^ + \v2\'^)} = Cl^/I. 

By Lemma 1 we also conclude that (44) holds and this completes the proof of Theorem 1. ■ 

4. Propagation of chaos 

In this section we extend the result obtained in Theorem 1 to the j-marginal distribution, 
showing the propagation of chaos (at first order in time). More precisely we have 

Theorem 2. Under hypotheses of Theorem 1, if (15) holds for all j, then 
(82) 



1=1 



J 

{t,xi,vi,. . .,Xj,Vj) = '[]_S{t)fo{xi,Vi)+ 

^'^ i=l 

^'[[S{t)fo{xk,Vk) / drSit - T)QLiS{T)fo,S{T)fo){xi,Vi), 

,—1 I 1 ■'O 



i=l k=l 

(83) lim'yj^ {t,Xi,vi, . . . ,Xj,Vj) = 
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in D'. 

Remark 2. The reason why we call eq.n (82) propagation of chaos is that the r.h.s of (82) 

corresponds to the first order in tim,e ofY\^^^^Ymig^Qgi[xi,iii,t). 



Proof. Let u G D(M'^-' x Mr^) be a test function and let us consider 

(u,fff(t)) = (n,5(i)/°) + — =^ / {u,S{t-T)q^,SiT)f^^,)dT 
(84) ^° 

+ ^ / (u,5(t-r)C-+i7f+i(r))dr. 

Of course the second term in (84) is of order 0{\/e), hence we focus on the third term. Then, 
defining Ri{Xj, Vj,T) := V^.5'(r — t)u{Xj, Vj), such a term is 



N-j 



jzfdrf / dy,+ii?,(X,,F,,r).F(^i^^)7f+i(r) 

Jo J J 



N-j 



= -^izf/r J dX,+, J dV,^, £ dsR,{X„V„r) g ^ (^^IZ^^ 



k<l 



{F (^^hLA_3LA^ . [(v., - V.J5(r-.)/jVi](X,+i(-.),lG+i(-.))}. 

We shall see that the leading term in the sum appearing in the r.h.s. of (85) is that with 
k = i,l = j + \, the other ones being vanishing. This is the content of the following 



Lemma 5. Let (p = (p{Xj j^i, Vj_|_i,t, s) > 6e a measurable function, compactly supported in 
Xj+i and such that 

Then, if {k, I) ^ {i,j + 1), for all i, k, I, we have 

(86) ^[dsj dX,,,jdV,^MF (^Z^) I IF (^^(Z£)^^) I < c^,. 

Proof. We are integrating on the final coordinates (Xj+i, V^+i) = (Xj+i(0), Vj-|_i(0)) of the 
flow (Xj_|_i((T), 1^_|_i((t)) defined for negative times a G [— r, 0]. We find convenient to reverse 
the velocities Vj+i — >■ — V^+i and look at positive times s G [0, r]. 

First of all we perform the usual change of variables a^j+i = xi — er and gain e^. Next 
we introduce the following partition of the phase space: setting Co = {(fc, Z),A; < l\{k,l) ^ 
+ 1)} we define 

(87) Aoik,l) = {iXj+i,Vj+i)\ \xk-xi\<2e,ik,l) gCo} 
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and 

(88) ^0= U Ao{k,l). 

Furthermore, denoting by s{k,l) G [0,r] the first instant for which 

(89) \xk{s) - xi{s)\ < e 

namely the pair of particles k and I starts to interact at time s{k, I) (if they do not interact 
we set s{k, I) = r) we define: 

(90) Ak,i = {{Xj+i,Vj+i) i ^o|s(fc, I) = , min s(r, m) < r}. 

(r,m)eCo 

In other words if (Xj+i, \^+i) G A^^i the pair of particles (A;,/) G Co is the first interacting 
pair (excluded the pair (i, j + 1) which starts to interact at time 0) in the time interval (0,r]. 
Note that we are interested to integrate over the set 

(91) AoU U ^M- 

(fc,OeCo 

In facts in the complement of the set (91), (86) vanishes because 

To estimate the contribution due to A^,; we first assume that A;/i,//j + l,i. 
Note that the motion of the pair of particles with indices (A;, Z) is free in [0, s(A;,Z)]. Then 
setting Xk — xi = y and — vi = w we have 

(92) inf \y — ws\ < e. 

se[o,T] 

The minimizing s is sq = so that condition (92) yields 

(93) |y_l(;!^m| <£. 

This means that the projection of y on the orthogonal plane to w is in the disk smaller than 
£. Therefore 

(94) ^ / dX,^,dV,^MF { ^i^^) I IF ( I < 

^ JA{k,l) V £ / V £ / 

Now we consider the cases k = i,l = ioYl=j + l. For the sake of clearness we consider 
k = i, the other cases being completely analogous. 

There are two possibilities: either s(i, I) > s, where s is the last interaction time for the 
pair + 1), namely 

\xi{s) - Xj+i{s)\ > e 

for s > s, or s(i, /) < s. 

In the first case we can repeat the above argument setting y = Xi{s) — xi{s) and w = 
Vi{s) - vi{s). 
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In the second one observe that the center of mass x = '^'+^-'+^ is moving freely with velocity 
V = ^'"'"^■'+^ (because the pair + 1) is an isolated system at least up to a time s = s{i, I)). 
Condition 

\Xi{s) - Xl{s)\ = £ 

implies 

3 

(95) \xi{s) - x{s)\ < \xi{s) - xi{s)\ + \xi{s) - x{s)\ < -e . 
Therefore we can integrate under the condition (95) to get 

(96) ^ / dX,^,dV,^, ^ \F ( ^l^^^l ||F ( Xk{-s)-x,{-s) \ I ^ ^^.^ _ 

^ JA{i,l) \ ^ J \ £ / 

Clearly we also have that 

(97) ^ 1^ dX,^,dV,^, ^ \F (^^izmi^ (^MZ£)_:^^ I < . 

Thus we conclude the proof. ■ 

Finally we handle the leading term. Setting 

rt 



Tl = -^^j^Y.J^ dT J dXj+i J dVj+1 £ dsRi{Xj,Vj,T)F 



(98) 
we have 

Lemma 6. The term with repeated indices is of order one. More precisely, 

( , 



(99) lim Ti 



3 J ,-t 
^^^^S[t)h{X},,Vk) \ dTS{t - T)QLiS{T)fQ,S{T)fQ){Xi,Vi),U 

1=1 k=l •'^ 



Proof. At this point the proof is rather obvious and we only sketch it. We first reduce the 
integration domain in the definition of Ti for moderately large relative velocity, i.e. \vi~Vjj^i \ > 
ae^/^, being the contribution of the complementary set negligible as we have seen in Section 
3. Looking at 

(100) 

i rt 

'-1 



£ 



I dT f dXj+i [ dVj+i r dsRi{Xj,Vj,ry 

Jo J J\vi-Vj+i\>ae'^/^ Jo 



we could apply the same argument as in Section 3 to get the result, if the motion of the pair 
of particles i and j + 1 would be independent of the others. However we have seen in the proof 
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of Lemma 5 that the contribution of the event in which the particle k ^ i,j + 1 interacts with 
particle i or particle j + 1 is indeed negligible. Hence (99) follows easily. ■ 

Finally, again by Lemma 1, we obtain (83). ■ 
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